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We present a non-thermal leptogenesis scenario following standard supersymmetric hybrid inflation, 
in the case where light neutrinos acquire mass via the usual seesaw mechanism and inflaton decay to 
heavy right-handed neutrino superfields is kinematically disallowed, or the right-handed neutrinos 
which can be decay products of the inflaton are unable to generate sufficient baryon asymmetry 
via their subsequent decay. The primordial lepton asymmetry is generated through the decay of 
the inflaton into light particles by the interference of one-loop diagrams with exchange of different 
right-handed neutrinos. The mechanism requires superpotential couplings explicitly violating a 
U(l) R-symmetry and R-parity. We take into account the constraints from neutrino masses and 
mixing and the preservation of the primordial asymmetry. We consider two models, one without 
and one with SU(2 )_r gauge symmetry. We show that the former is viable, whereas the latter is 
ruled out. Although the broken R-parity need not have currently observable low-energy signatures, 
some R-parity-violating slepton decays may be detectable in the future colliders. 
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I. INTRODUCTION 

One of the most promising scenarios for generating the 
observed baryon asymmetry of the universe (BAU) is cer¬ 
tainly the leptogenesis scenario Hi It applies in all the 
cases where the light neutrinos {y) acquire their mass by 
coupling to heavy standard model (SM) singlet fermions 
v'^, the right-handed neutrinos (RHNs) (this is known 
as the seesaw mechanism Q)> or SU(2)i triplet Higgs 
scalars 0. These heavy particles can decay out of ther¬ 
mal equilibrium generating a primordial lepton asymme¬ 
try, which is subsequently converted in part into the ob¬ 
served baryon asymmetry by non-perturbative sphaleron 
effects at the electroweak phase transition. 

In the original realization H of this scenario, the heavy 
particles were assumed to be thermally produced in the 
early universe. However, there is a tension between cor¬ 
rect neutrino masses and this thermal leptogenesis sce¬ 
nario in supersymmetric (SUSY) models because of the 
gravitino problem mi. Assuming that the gravitino 
mass is of the order of 1 TeV and employing generic 
sparticle spectra, one finds that the reheat temperature 
Treh should not exceed about 10® GeV since otherwise 
an unacceptably large number density of gravitinos is 
thermally produced at reheating. These gravitinos later 
decay presumably into photons and photinos interfering 
with the successful predictions of standard big bang nu¬ 
cleosynthesis. On the other hand, adequate thermal pro¬ 
duction of RHNs or SU(2 )l triplets, whose decay creates 
the primordial lepton asymmetry, requires that the mass 
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of these particles does not exceed Treh. This leads to un¬ 
acceptably large light neutrino masses. This problem can 
be alleviated H @1 if we allow some degree of degeneracy 
between the relevant RHNs, which enhances the gener¬ 
ated lepton asymmetry, and perhaps also if the branching 
ratio of the gravitino decay into photons and photinos is 
less than unity, which somewhat relaxes the gravitino 
constraint on Treh- 

A much more natural solution of the tension between 
the gravitino bound on Tj-eh and the masses of light neu¬ 
trinos is provided by non-thermal leptogenesis Q at re¬ 
heating. In existing realizations of this scenario, 
though, where the inflaton decays into RHN or SU(2)j;, 
triplet superfields, this still puts a restriction on the 
masses of these particles: the decay products of the in¬ 
flaton must be lighter than half its mass minf. The pri¬ 
mordial lepton asymmetry is generated in the subsequent 
decay of the RHN or SU(2)i triplet superfields. 

In a recent paper El, we considered the consequences 
of allowing all the RHN and SU(2 )l triplet superfields 
to be heavier than minf/2 (see also Ref. 0 ). Primor¬ 
dial leptogenesis could then take place only through the 
direct decay of the inflaton into light particles (see also 
Ref. ^11). We took a simple SUSY grand unified theory 
(GUT) model which is based on the gauge group Gb-l = 
Gsm X U(l) b-l (Gsm is the SM gauge group, and B and 
L the baryon and lepton number respectively) and natu¬ 
rally incorporates the standard SUSY realization Ellll 
of hybrid inflation M- The flatness of the inflation¬ 
ary trajectory at tree level was guaranteed by a U(l) R- 
symmetry, whereas radiative corrections provided a 
logarithmic slope along this path, needed for driving the 
inflaton towards the SUSY vacua. The R-symmetry also 
guaranteed the conservation of baryon number to all or¬ 
ders in perturbation theory. Therefore, baryon number 
was only violated by the non-perturbative electroweak 
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sphaleron effects. 

The model incorporated the solution of the ^ problem 
of the minimal su per symmetric standard model (MSSM) 
proposed in Ref. [^. Although the global R-symmetry 
forbade the appearance of a /r term in the superpoten¬ 
tial, it did allow the existence of the trilinear term Shih 2 , 
where S is the gauge singlet inflaton of standard SUSY 
hybrid inflation and hi, /12 are the electroweak Higgs su¬ 
perfields. After the GUT gauge symmetry breaking, the 
soft SUSY-breaking terms, which generally violated the 
R-symmetry, gave rise to a suppressed linear term in S 
and, thus, this field acquired a vacuum expectation value 
(VEV) of the order of the electroweak scale divided by 
a small coupling constant. The above trilinear coupling 
could then yield a /i term of the right magnitude. 

The same coupling also gave rise to tree-level couplings 
of the inflaton, which consisted of two complex scalar 
fields, to the electroweak Higgs bosons and Higgsinos. 
After the termination of inflation, the inflaton performs 
damped oscillations about the SUSY vacuum and even¬ 
tually decays predominantly into electroweak Higgs su¬ 
perfields via these tree-level couplings, thereby reheating 
the universe. The model contained both heavy RHN and 
SU(2 )L triplet superhelds which both contributed to light 
neutrino masses. A primordial lepton asymmetry could 
be generated at reheating via the subdominant decay of 
the inflaton into lepton and Higgs superhelds through the 
interference between one-loop diagrams with RHN and 
SU(2)i triplet exchange respectively. The simultaneous 
presence of both RHNs and SU(2 )l triplets was essen¬ 
tial for this particular leptogenesis mechanism to work. 
Note that, in the model of Ref. El. the generation of a 
non-zero lepton asymmetry did not rely on the existence 
of more than one fermion family. However, more than 
one family was required if the lepton asymmetry was to 
be preserved down to the electroweak phase transition. 
Finally, it should be emphasized that the generation of 
a non-zero lepton asymmetry required the inclusion of 
some couplings in the superpotential that explicitly vio¬ 
late the U(l) R-symmetry. 

In this paper, we investigate the consequences of doing 
without SU(2 )l triplet superhelds and instead generating 
the primordial lepton asymmetry through the exchange 
of different RHN superhelds. Since these heavy helds can 
only appear in intermediate states of the inhaton decay, 
we must create the asymmetry directly from this decay. 
Indeed, leptogenesis can again occur in the subdominant 
decay of the inhaton into lepton and Higgs superhelds, 
but now through the interference between one-loop dia¬ 
grams with diherent RHN exchange. The lepton asym¬ 
metry is proportional to a novel CP-violating invariant 
product of coupling constants. 

We study this new leptogenesis scenario within the 
framework of the model of Ref. 0 with the SU(2 )l 
triplet superhelds removed. All the salient properties of 
the model are retained except that now the light neu¬ 
trino masses are generated by the standard seesaw mech¬ 
anism 0 which involves only RHN (super)helds and lep¬ 


togenesis takes place via the new mechanism mentioned 
above. In particular, the implementation of the standard 
SUSY hybrid inhationary scenario MM and the so¬ 
lution of the /i problem of Ref. 0 remains unaffected. 
Also, baryon number is still conserved to all orders in 
perturbation theory. 

The model again contains an approximate U(l) R- 
symmetry explicitly broken by some superpotential op¬ 
erators. These operators, which are necessary to create 
a non-zero lepton asymmetry, also violate the Z 2 matter 
parity subgroup of the R-symmetry which is preserved 
by soft SUSY-breaking terms. The matter parity viola¬ 
tion may have important observable consequences at low 
energy. Indeed, if the lightest sparticle (LSP), which is 
assumed to be the lightest neutralino, contains a Hig- 
gsino component, we hnd that it could decay predomi¬ 
nantly into one or two Higgs bosons and a lepton. These 
dominant decay channels can though be easily blocked 
kinematically if the LSP is not too heavy and, under cer¬ 
tain conditions, this particle can be made long-lived (see 
Sec.IVn. However, some matter-parity-violating slepton 
decays which may be detectable in the future colliders 
are typically present (see Sec. IIXII . 

We perform a detailed chemical potential analysis of 
the evolution of the primordial lepton asymmetry un¬ 
til the time at which the baryon-number-violating elec¬ 
troweak sphaleron effects cease to operate and evalu¬ 
ate the baryon asymmetry at this moment, which yields 
the observed BAU. We find that there are some lepton- 
number-violating four-scalar processes involving Higgs 
bosons and sleptons which are in equilibrium in a range 
of temperatures just above the SUSY threshold. These 
processes (which result from the same dimension-four 
operators that cause the above-mentioned slepton de¬ 
cays) could readily erase the primordial lepton asym¬ 
metry, which would lead to the absence of any baryon 
asymmetry in the present universe. However, it is possi¬ 
ble to choose parameters so that one of the three lepton 
family numbers is conserved by all dimension-four pro¬ 
cesses. One can further show that all the dimension-five 
processes which violate this number are well out of equi¬ 
librium at all temperatures after reheating. Thus, we can 
obtain a non-zero BAU at present. 

We find that the value of the BAU from the Wilkinson 
microwave anisotropy probe (WMAP) data can be 
achieved given constraints from other observables, no¬ 
tably the reheat temperature and neutrino masses and 
mixing, and CP-violating phases of order unity. Note 
that, contrary to Ref. where we worked in the two 
heaviest family approximation, in the present paper we 
take fully into account all three neutrino species. We 
find that only for a small range of values of the lightest 
neutrino mass mi can the scenario be successful, given 
the requirement that the GUT breaking scale be pertur¬ 
bative relative to an underlying string or quantum grav¬ 
ity scale, which is also restricted to smaller values than 
the value used in Ref. 0 . The reheat temperature is 
pushed to higher values relative to this reference. Finally, 
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the values of the neutrino Yukawa couplings are also re¬ 
stricted. These results reflect the fact that removing the 
SU(2 )L triplets from the model makes it more restrictive. 
The prediction for the spectral index of density pertur¬ 
bations is typical of SUSY hybrid inflation models (see 
e.g. Ref. O). 

Thus, an acceptable value of the baryon asymmetry 
can be obtained within a consistent model of cosmology 
and particle physics, without requiring additional fine- 
tuned coupling constants, while respecting experimen¬ 
tal constraints on neutrino masses and mixing. More¬ 
over, although the scenario requires violation of the R- 
symmetry, it is not necessary to introduce superpoten¬ 
tial terms which would lead to currently observable R- 
symmetry-violating effects. 

We also analyzed the case where SU(2)i{ gauge symme¬ 
try is imposed, forcing the neutrino Yukawa couplings to 
be identical to those of the charged leptons at the GUT 
scale. This case is thus considerably more restrictive than 
the non-SU(2)/i-symmetric case. One of the main differ¬ 
ences is that it predicts that the lightest RHN mass is 
of order 10^ GeV, whereas the inflaton mass is greater 
than 10^^ GeV, thus the inflaton can decay directly to 
the lightest RHN. However, this does not alter our basic 
picture for leptogenesis since the corresponding branch¬ 
ing ratio and non-thermally generated lepton asymmetry 
are negligible due to the fact that the lightest RHN mass 
is very small compared to the inflaton mass. Moreover, 
thermal leptogenesis from the decay of this particle can 
be ignored for the same reason. The presence of the 
lightest RHN in the thermal bath after reheating, how¬ 
ever, yields extra restrictions on the parameters of the 
model from the requirement that the initial lepton asym¬ 
metry is not erased. In the end, this case is ruled out be¬ 
cause it always yields an unacceptably small BAU. The 
reason is that our leptogenesis scenario requires strong 
mixing in the RHN sector but, in this case, the RHN 
mass matrix has too small off-diagonal elements due to 
the strong hierarchy of Yukawa coupling constants im¬ 
plied by the SU(2)fl symmetry and the extra restrictions 
from the preservation of the primordial asymmetry. 

In Sec. m we introduce the SUSY GUT models and 
describe some of their salient features. In Sec. nni we 
discuss constraints from the light neutrino masses and 
mixing, while, in Sec. CYI we present the CP-violating 
invariant products of coupling constants which enter into 
the primordial lepton asymmetry. The calculation of the 
primordial lepton asymmetry is sketched in Sec. m and 
the effects of R-symmetry (and R-parity) violation are 
discussed in Sec. ED Section Em is devoted to the con¬ 
ditions for the primordial lepton asymmetry to be pre¬ 
served at high temperatures. We revisit the relation be¬ 
tween the initial lepton asymmetry and the final baryon 
asymmetry in Sec. IVHII finding some novel results for 
the case when some lepton family numbers are violated. 
Our numerical results appear in Sec. im including a dis¬ 
cussion of why the SU(2)/j-symmetric case is ruled out. 
Finally, conclusions appear in Sec. HI 


TABLE I: U(l) charges of superfields 



S 

0 


hi 

h2 

1 



q 



B-L 

0 

1 

-1 

0 

0 

-1 

1 

1 

1/3 

-1/3 

-1/3 

R 

2 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 


II. THE SUSY GUT MODELS 

We consider two SUSY GUT models which are based 
on the gauge groups Gb-l and Glr = SU(3)cXSU(2)lX 
SU(2)kxU(1)b_l respectively. The models also possess 
an approximate global R-symmetry U(l)/{, which is ex¬ 
plicitly broken by some terms in the superpotential. 

We will treat the Gb-l case first. In addition to the 
usual MSSM chiral superfields hi, h2 (Higgs SU(2)i dou¬ 
blets), li (SU(2)i doublet leptons), ef (SU(2 )l singlet 
charged leptons), qi (SU(2)i doublet quarks), and uf, dj 
(SU(2 )l singlet anti-quarks) with i = 1,2,3 being the 
family index, the model also contains the SM singlet chi¬ 
ral superfields (RHNs), S, (j), and (f>. The charges 
under \J{1)b-l and U(l)fl are given in Table |D The 
superpotential is 

W = M'^) + \S{hih2) + yeijihhi)ej 

+yuijiqth2)uj + ydrj{qihi)dhj 
+yi^ij{kh2)Rj + 

-\-{Xi/Ms)^i'1{hih2) + ■ ■ ■ , ( 1 ) 

where M is a mass parameter of order the GUT scale, Ms 
is the string or quantum gravity scale ^ 10^”^ GeV, and 
(XY) indicates the SU(2)i invariant product eabXaYh 
with e denoting the 2x2 antisymmetric matrix with 
ei 2 = 1. The ellipsis represents terms of order higher 
than four and summation over indices is implied. The 
only U(l)/{-violating couplings which we allow in the su¬ 
perpotential are the ones in the last explicitly displayed 
term in the right hand side (RHS) of Eq. (E), which are 
necessary for leptogenesis. It can be shown that baryon 
number (R) is automatically conserved to all orders as a 
consequence of U(l)/{. The argument goes as in Ref. M 
and is not affected by the presence of the above U(l)/{- 
breaking superpotential couplings. Lepton number (L) is 
then also conserved in the superpotential as implied by 
the presence of U(1 )b_l. 

In the SU(2)/{-symmetric case, the field content is sim¬ 
plified as follows: the SU(2)i singlets and form the 
SU(2)/{ doublets l^, the anti-quark superfields and 
form the SU(2)/{ doublets qf, and the electroweak 
Higgs superfields hi and h2 form the SU(2 )l x SU(2)/{ 
bi-doublet h = {h2,hi). The symmetry-breaking sin¬ 
glets (j) and ([> are replaced by SU(2)/{ doublets M and 
whose neutral (RHN-like) components N'^ and obtain 
VEVs of order the GUT scale. To simplify notation, we 
will denote these components also by (j) and (j) respec¬ 
tively. The superfields M and M have the same charges 
under \]{1)b-l and U(l)i{ as 4> and (j) respectively, while 
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the charges of ql and h coincide with the charges of 
their components in Table Q] Below the scale M, the 
field content is the same as in the previous case. The 
superpotential, in this case, is 

W = KS(L^L‘^-M^) + ASh^ + yHjhehl^ 

+yq^Jq^ehq‘^ + 

+ {K/Ms)L%h^+ ■■■, ( 2 ) 

where denotes the gauge invariant sum iTr(eft.^e/i) = 
(/ 11 / 12 ) and i'^, If, qf are taken as column 2 -vectors, while 
i'^, li, qi as row 2-vectors (the superscript T denotes ma¬ 
trix transposition). Again the only 11(1)^-violating su¬ 
perpotential couplings are the ones in the last explicitly 
displayed term in the RHS of Eq. and B conserva¬ 
tion to all orders is automatic. Note that, imposing the 
SU(2)i{ symmetry on the superpotential of Eq. m , we 
obtain that the Yukawa coupling constants yeij {yuij) and 
Vvij iudij) become equal |^. Their common value is the 
coupling constant yuj (jjqij) in Eq. Q. 

In both cases, the classically flat direction in field 
space along which inflation takes place is as described in 
Ref. for K < A, it is parametrized by S', [S'! > S'c = 
M, with the values of all the other fields being equal to 
zero, and has a constant potential energy density 
(at tree level). Here, the dimensionless coupling con¬ 
stants K, X and the mass parameter M are taken real and 
positive by appropriately redefining the phases of the su¬ 
perfields. There are radiative corrections |l fJl which lift 
the flatness of this classically flat direction leading to 
slow-roll inflation which, for k 1 , terminates practi¬ 
cally when |S| reaches the instability point at IIS'! = M 
as one can easily deduce from the slow-roll e and 77 cri¬ 
teria M The quadrupole anisotropy of the cosmic mi¬ 
crowave background radiation (CMBR) and the number 
of e-foldings 


Nq ~ In 


1.88 X 




(3) 


suffered by our present horizon scale during inflation are 
given by the Eqs. (2)-(4) of Ref. |^; in the non-SU(2)fl;- 
symmetric case, the two last terms in the RHS of Eq. (3) 
must be divided by two since the SU(2)fl doublet chiral 
superfields are replaced by the SM singlets (j), ^ 

(see also Ref. I^h 

For K < A, a GUT-symmetry-breaking phase transi¬ 
tion takes place when the value of [S'! falls below the 
mass parameter M. The fields evolve towards the re¬ 
alistic SUSY minimum at (S) = 0, {(j)) = {^) = M, 
{hi) = (/ 12 ) = 0, where the GUT-symmetry-breaking 
A^EVs are taken real and positive by a i? — L rotation. 
There is also an unrealistic SUSY minimum which is 
given below. With the addition of soft SUSY-breaking 
terms, as required in a realistic model, the position of 
the vacuum shifts to non-zero (S) ~ where 

m 3/2 is the mass of the gravitino, and an effective y term 


with n ~ —Xm^/ 2 lK, ~ m ^/2 is generated from the super¬ 
potential coupling XS{hih 2 ). 

After this phase transition, the (complex) inflaton de¬ 
grees of freedom are S and 9 = {S(j) + 5(j))/i/2, where 
d(j) = (j) — M and 5(j) = (j) — M, with mass minf = V^kM. 
These fields oscillate about the minimum of the poten¬ 
tial and decay to MSSM degrees of freedom reheating 
the universe. The predominant decay channels of S* 
and 9 are to fermionic and bosonic hi, /12 respectively 
via tree-level couplings derived from the superpotential 
terms XS{hih 2 ) and kS4>4> (or kSL‘^L‘^). Note that, if 
At > A, the system would end up in the unrealistic SUSY 
minimum at </) = 0 = 0, |/ii| = II 12 I — {k/X)^^‘^M, which 
is degenerate with the realistic one (up to 171 ^ 2 ) is 
separated from it by a potential barrier of order 

The RHNs acquire masses after the spontaneous 

breaking of U(1 )b_l by {(j)), {(j)). The coupling constants 
of the relevant terms in the third line of the RHS of 
Eq. (or the second line of Eq. iPJ) can also be writ¬ 
ten as Xv^ij/Ms, making it clear that the RHN masses 
are suppressed by a factor M/Ms relative to M. It is 
possible to redefine the superfields I'f (or If) to obtain 
effective mass terms Mivfvf which are diagonal in the 
flavor space with Mi real and positive. This basis, which 
we will call “RHN basis”, is most convenient for calcu¬ 
lating the BAU (see Sec. 0. For definiteness, we take 
Ml < M 2 < M 3 . 

We will also use a basis in which the lepton Yukawa 
couplings and the SU(2 )b interactions are diagonal in fla¬ 
vor space, denoting couplings and fields in this “lepton 
family basis” with a hat. We thus have yeij = Sijjjei, 
yuij = Sijiji (or ijiij = dijitii). The diagonal Yukawa 
coupling constants yei, IJi (or yu) are taken real and pos¬ 
itive by appropriate rephasing of the fields. Note that, 
in the non-SU(2)B-symmetric case, it is an additional as¬ 
sumption that the neutrino and charged lepton Yukawa 
couplings can be diagonalized in the same weak inter¬ 
action basis. We will elaborate on this assumption in 
Sec, ivm Since the weak interactions are in equilibrium 
at temperatures above the critical temperature Tc for the 
electroweak phase transition, the lepton family numbers 
Li can only be defined in this basis. 

Relative to the “unhatted” RHN basis, we have 

i=lU, = (4) 

Here U, C/“ are 3x3 unitary matrices and we write 
“left-handed” lepton superfields, i.e. SU(2)b doublet lep¬ 
tons, as row 3-vectors in family space and “right-handed” 
anti-lepton superfields, i.e. SU(2)b singlet anti-leptons, 
as column 3-vectors. In the SU(2)B-symmetric case, we 
have the same relation except that the SU(2 )b doublet 
B (L) replaces iy‘^ (i^'^)- Then by definition we have, for 
the neutrino Yukawa couplings, 

leh2y^i'‘' = ieh2U^y^U‘'^i>^, 

U^yuU’'^ = = diag(yi,y2,y3), 


(5) 
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where yu is the 3x3 matrix with elements yuij- In 
general, the jji can be any (real positive) numbers. In 
the SU(2)ij-symmetric case, however, they are deter¬ 
mined by the “asymptotic” values (at the GUT scale) of 
the charged lepton masses as j/ 1 , 2,3 = We,/j,r tan/3/(/i2)- 
The Higgs VEV is (/ 12 ) — 174sin/3 GeV, which, in the 
large tan/3 limit, yields (/ 12 ) ~ 174 GeV. Exact SU(2)i{ 
Yukawa coupling relations in the quark sector would fix 
the ratio of the Higgs VEVs, but at least small deviations 
are required here to be consistent with the data. Assum¬ 
ing that the corrections are small compared to the third 
generation Yukawa couplings, tan/3 should be about 55. 
(For the conditions under which SU(2)/j implies large 
tan/3, see Ref. [l^.) The matrices U and C/'^ are, at this 
stage, determined only up to a diagonal matrix of ar¬ 
bitrary complex phases P = diag(e*‘^b 6*“^^) which 
acts as U ^ UP, U‘^ P~^U‘^, corresponding to oppo¬ 

site phase redefinitions of the lepton weak doublet and 
singlet superfields. 

As for the U(1 )_r- violating terms, after the U(1 )_b_i, 
(and SU(2)/{) breaking, the explicitly displayed terms 
in the last line of the RHS of Eq. o or 0 give rise 
to effective B — L and matter-parity-violating operators 
Cit'f(hi/ 12 ), where the dimensionless effective coupling 
constant Q is suppressed by one power of M/Ms- If 
we require that the magnitude of the dimensionless cou¬ 
pling constants Xi be less than unity, we obtain a bound 
ICi| ^ M/Ms- Note that the effective coupling constants 
Ci cannot in general be made real and positive by redefin¬ 
ing the complex phases of the superfields. This can be 
easily shown by considering the phase redefinition invari¬ 
ant y?Mi (no summation over repeated indices) with 
p. and Mi already made real (recall that the Q are in the 
RHN mass eigenstate basis). The coupling constants Q, 
thus, remain in general complex. Note that {hi), (/ 12 ) can 

many other 

R-symmetry-violating operators. However, they are not 
necessary for a non-zero primordial lepton asymmetry to 
be created. 

The Yukawa coupling constants yuij (in the original 
RHN mass eigenstate basis) remain in general complex 
as one can easily deduce from the rephasing invariants 
VvijVtikCjCk for j ^ k (no summation over repeated 
indices), as well as the standard rephasing invariant 
yuijyukjytijnVtkmMjMm for j ^ m (no summation), 
which has been used in the original leptogenesis scenario 
of Ref. 0. 

The calculation of lepton asymmetry produced in S 
and 9 decays is quite straightforward but somewhat 
lengthy, and differs in detail from the usual case where 
leptogenesis takes place through the decay of RHN or 
SU(2)i triplet superfields. Since we consider the inter¬ 
ference of two one-loop diagrams, we will need to find 
the real parts of loop integrals, which require renormal¬ 
ization. Before proceeding to the calculation, we sum¬ 
marize the constraints derived from the current neutrino 
experiments. 


be taken real because of the reality of By ~ 
Ea. It is, of course, possible to write down 


III. NEUTRINO MASSES AND CONSTRAINTS 

The standard seesaw mechanism yields the light neutrino 
mass matrix: 


m = - {h 2 )'^yuM,^}yl, (6) 

which holds in any basis of neutrino states with 
being the 3x3 matrix with elements M^^ij. In the RHN 
mass eigenstate basis, where = diag(Mi, M 2 , M 3 ), 
it can be rewritten as 


m = -{h2)^Uy,U^MU.^U^^y.U^. (7) 

Thus, transforming the light neutrino fields to the “hat¬ 
ted” (lepton family) weak interaction basis, we have 

m = UlmU* = -{h2)^yuU‘'Mu}U‘^'^y,,. (8) 

The light neutrino mass eigenstate basis can be denoted 
by V and is given by 


D = vV*, (9) 

where v, v are row 3-vectors in family space and V is 
the unitary Maki-Nakagawa-Sakata (MNS) matrix which 
satisfies 


= m = diag(TOi, m 2 , TO 3 ) ( 10 ) 


with rrii real and positive. Again V is determined only 
up to a diagonal matrix of complex phases as U ^ PV, 
which reflects the freedom to redefine the phases of the 
superfields in the “hatted” basis. However, the ambiguity 
is fixed by requiring V to take the standard form 


V = 


C 3 C 2 

-C1S3 — C3SiS2e 


iS 


S 3 C 2 
C3C1 — S3SiS2e 


— i6 


iS 


S2e 

C2S1 


S3S1 - C3CiS2e*'^ -C3S1 - CiS3S2e*^ C2C1 

( 11 ) 

where Ci = cos 9i, Si = sin^i with 9i {i = 1,2,3) being 
the appropriate mixing angles, S is the Dirac phase and 
V = diag(e“®“, 1) contains the Majorana phases. 

In both the SU(2)/j- and non-SU(2)/j-symmetric cases, 
we will restrict the light neutrino mass-squared differ¬ 
ences and mixing angles to be within the 2 a allowed re¬ 
gion detailed in Ref. j^. Eor most numerical results, 
however, we will set these observables to their best-fit val¬ 
ues, implying that m is completely determined once the 
lightest neutrino mass, the phases a and (3 and the choice 
of normal or inverted hierarchy are fixed. We will take 
here the light neutrino mass ordering mi < m 2 < m 3 
and, unless stated otherwise, adopt the normal hierarchi¬ 
cal scheme of neutrino masses, where the solar and atmo¬ 
spheric neutrino mass-squared differences are identified 
with (5m|i and (5m§i respectively (here 5m1j = mi—rrij). 
For simplicity, we further take a = /3 = 0. 

Eor given values of m^ and V, one can find the light 
neutrino mass matrix in the “hatted” basis by invert¬ 
ing Eq. cni). Then, in the SU(2)/{-symmetric case, the 
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seesaw formula in Eq. m can be applied with iji deter¬ 
mined by the known values of the charged lepton masses, 
renormalization group (RG) evolved to the GUT scale, 
and for given tan/3. Solving the resulting equations, we 
can simultaneously determine the elements of the unitary 
matrix U‘^ and the RHN mass eigenvalues Mi. In the non- 
SU(2)/{-symmetric case, the Yukawa coupling constants 
yi are not related to the charged lepton masses and, thus, 
are treated as free input parameters. 

For our mechanism to be the exclusive source of baryon 
asymmetry as desired, either all RHN masses must be 
greater than minf/2 so that the inflaton decay to RHN 
superfields is kinematically blocked and the thermal pro¬ 
duction of RHNs after reheating is suppressed, or any 
RHNs that are lighter must contribute a negligible net 
lepton number density via the conventional thermal or 
non-thermal leptogenesis mechanism, by which they de¬ 
cay out of equilibrium. One way to achieve the latter is 
to make them adequately light. We also have to check 
that thermal interactions of any lighter RHNs do not 
wash out the already created lepton asymmetry. Fur¬ 
thermore, assuming that the dimensionless coefficients 
of the last superpotential terms in the third line 
of the RHS of Eq. (or the second line of the RHS of 
Eq. (|21)) do not exceed unity, we obtain that the heaviest 
RHN must not have a mass greater than /Ms- These 
requirements place non-trivial constraints on the model 
parameters which will be further discussed in Sec. nsi 


IV. CP-VIOLATING INVARIANTS 

The generation oi a B — L asymmetry from the decay 
of the inflaton requires that the theory contains at least 
one CP-violating product of coupling constants which re¬ 
mains unaffected by field rephasing and is non-real cor¬ 
responding to an operator non-invariant under CP con¬ 
jugation. Since leptogenesis takes place at reheating, we 
work in the vacuum where ((jX) = ((/>) = M and ignore 
the /i term. In previous work |ll|. when calculating the 
lepton asymmetry from the CP-violating inflaton decay, 
we considered only the couplings of the inflaton originat¬ 
ing from the renormalizable terms in the superpotential 
of Eq. (or Eq. ||2Il). In fact, when one does also in¬ 
clude the couplings of the inflaton to RHN and Higgs 
superfields from the last two explicitly displayed non- 
renormalizable terms in the RHS of Eq. o (or Eq. @), 
it turns out that their contribution to the asymmetry 
vanishes exactly. Thus, we will replace these two super¬ 
potential terms by the effective mass terms and 

couplings Ci’^f{hih 2 ) respectively. Since the only massive 
particles at reheating are the inflaton and the RHNs, it is 
convenient to work in the basis where the RHN Majorana 
masses are diagonal. 

In Ref. 0 , we considered rephasing invariants which 
could survive even if there were only one matter gener¬ 
ation. This was possible because of the presence of the 
SU(2)/:, triplet superfields. In the current models, we 


must rely on the presence of more than one generation. 
In addition to the standard leptogenesis CP-violating in¬ 
variant composed of four Yukawa coupling constants and 
two RHN masses, there is a CP-violating invariant which 
can be written in the “hatted” basis, where the neutrino 
Yukawa couplings are diagonal, as 

lo^jk = ( 12 ) 

for i ^ j (no summation). However, this does not survive 
transformation to the RHN basis and does not on its own 
give rise to diagrams which can create a lepton number 
asymmetry in light fields since the corresponding opera¬ 
tor does not involve any light lepton superfields. Instead, 
we “dress” it with Yukawa coupling constants to obtain 
the CP-violating invariant 

hjk = (13) 

for j ^ k (no summation), where we now use the mass 
eigenstate basis of RHNs. When this invariant is used, 
the generation oi a B — L asymmetry is independent of 
the sources of CP violation considered in previous sce¬ 
narios and we require novel decay diagrams. This CP- 
violating invariant is minimal in the sense that it involves 
the smallest possible number of trilinear superpotential 
couplings. We can transform the first (light lepton) index 
of this invariant to the “hatted” basis and re-express the 
invariant in terms of the “hatted” coupling constants as 

4fc = E (m) 

m,n 

where (m is defined by 

Cn=E4c^™, (15) 

m 

and we have redefined fields such that Mi and yi are real 
and positive (the two last indices of the invariant remain 
in the RHN basis). This will be convenient because the 
Yukawa couplings in this basis preserve the individual 
lepton numbers Li. 

Note that the invariant Ljk can be split in two factors 
yvijCjMj and ytikCk^^k which correspond to effective op¬ 
erators with opposite non-zero B — L charges. These 
operators involve only light fields since the heavy ones 
can be contracted. One of these two operators includes 
bosonic or fermionic hi, h 2 , while the other contains their 
conjugates. These properties are essential for leptogene¬ 
sis since the inflaton field couples at tree level to the elec- 
troweak Higgs superfields hi, h 2 and can decay only to 
light particles. So, we conclude that Ljk is, in principle, 
suitable for leptogenesis which requires the interference 
of two {B — L)-violating diagrams for the inflaton decay. 
One can further show that, if the RHN spectrum and 
couplings are such that the standard leptogenesis sce¬ 
nario (without the C, couplings) fails, any invariant which 
can be useful for leptogenesis must involve lijk, and thus 
the effective coupling constants Q. So, the explicit vi¬ 
olation of U(l)/{ and matter parity is essential for our 
scheme. This is another novel feature of this scenario. 
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FIG. 1: The six one-loop diagrams for the Li-violating decay 
S —> lih 2 . The solid (dashed) lines represent the fermionic 
(bosonic) component of the indicated superfield. The arrows 
depict the chirality of the superfields and the crosses are mass 
insertions in fermion lines. 


V. BARYON ASYMMETRY 

In Figs. Hand 12 we present the diagrams for the Li- 
violating decay of S and 9* respectively. Observe that 
the CP-violating rephasing invariant lijk corresponds to 
the product of coupling constants in the interference of 
any diagram in Fig. ^ (Fig- lit with any diagram in the 
same figure where the exchanged RHN is renamed 
This interference contributes to the Li asymmetry due 
to a partial rate difference in the decays S ^ lih 2 and 
S* I* /i 2 {0* lih 2 and 9 lih 2 ), where bar and 
tilde represent the anti-fermion and the SUSY partner 
respectively [2^. Both the decaying inflaton field {S or 
9*), which is taken at rest, and the decay products must 
be on mass shell. For simplicity, we consider that all 
the propagating and external MSSM particles in the di¬ 
agrams are massless. Moreover, we perform the calcula¬ 
tion in the limit of exact SUSY. 

We also considered diagrams where the inflaton de¬ 
cays via its couplings to RHN and Higgs superfields 
which arise from the superpotential terms with coeffi¬ 
cients Mi,cij and A^. Such diagrams exist only for the 

decay of 9* and have final states {h hi, k hi, k hi, I* hi) 
which are different from the final state of the diagrams 
shown in Fig. 12 However, it can be argued that these 


FIG. 2: The two one-loop diagrams for the Li-violating decay 
9* —> U h 2 . The notation is as in Fig.Q 

diagrams sum to zero in the decay amplitudes if the fi¬ 
nal state is khl or I* hi. For the other two cases, the 
diagrams again sum to zero but only in the amplitudes- 
squared after summing over the phase space of final par¬ 
ticles and in the limit where these particles are massless. 

We will denote by and the “stripped” di¬ 

agrams in Figs. and^ respectively with the dimen¬ 
sionless coupling constants and the Mj mass insertions 
factored out (we keep, though, the mini factor appearing 
in the scalar coupling 9 *hih2). Here i and j are family 
indices and n = 1,2,3 the serial number of the diagram. 
Similarly, the “stripped” diagrams in Figs.[2tt and|23 are 
F^ij and F^^j. In each case, the contribution to the Li 
asymmetry is proportional to both Im lijk and the imag¬ 
inary part of the interference of the relevant “stripped” 
diagrams. Thus, the total net Li asymmetries ei\s and 
generated per S and 9* decay respectively are 

j,k ,n,n' 

uie = ( 16 ) 

j.,k t.,t' 

where F = |ApTOinf/87r is the common rate of the tree- 
level decays S ^ hi 9,2 and 9 ^ hi / 12 , the indices t, t' = 
a, b, and integration over the phase space of the particles 
in the final state is implied. 






The diagrams are both of the vertex (Figs.^ andl^t-) 
and self-energy (Figs.^^ and|2t>) type. Each of the 
three vertex diagrams in Fig. QJi possesses a logarithmic 
ultraviolet (UV) divergence. However, it can be easily 
shown that their sum equals minf times the vertex dia¬ 
gram in Fig. Hl, which is UV finite. Similarly, one can 
show that the sum of the three quadratically divergent 
self-energy diagrams in Fig.^ equals minf times the self¬ 
energy diagram in Fig. Wp The latter is though not UV 
finite. It rather possesses a logarithmic divergence and, 
thus, needs renormalization. 

The above relations between the diagrams for the Li- 
violating decays of S and 9* imply that e^s = ^i\e = ^i- 
We thus concentrate on the calculation of e^g. The ver¬ 
tex diagram in Fig. being finite (both its real and 
imaginary parts) is independent of the renormalization 
scheme used. However, the diagram in Fig.[2lD involving 
a divergent self-energy loop requires us to apply a renor¬ 
malization condition. As argued in Ref. 0 , the appropri¬ 
ate renormalization scheme, in this case, is the on-shell 
(OS) scheme. In a general theory with scalars Si, the OS 
conditions on the renormalized self-energies Ily(p^) are 
as follows: 


Re Hy (fij) = Re Hy (^|) = 0 (17) 

for the off-diagonal self-energies {i ^ j) and 

lim — 2 R-enii(p^) = 0 (18) 

for the diagonal ones (see e.g. Ref. US])- Here we take 
a basis where the renormalized mass matrix is diago¬ 
nal with eigenvalues fj,i. The imaginary part of the self¬ 
energies is finite and, thus, not renormalized. In Fig. |21 d, 
we have an off-diagonal self-energy diagram between the 
scalars 6 and i>|. Given that 9 is on mass shell, the real 
part of this diagram vanishes in the OS scheme. The 
imaginary part, however, gives a finite contribution. 

The asymmetry generated in the lepton number Li per 
decaying inflaton is then 


This equation for the contribution of the self-energy di¬ 
agram holds 0 provided that the decay width of is 
l^?nf ~ which is well satisfied in our model 

if Mj is not unnaturally close to minf. Also, note that the 
contribution to the second sum in the RHS of Eq. ra 
originating from the interference of two self-energy dia¬ 
grams (ie. the term with t = t' = self) vanishes. The 
reason is that is pure imaginary. Finally, the con¬ 
tributions to the RHS of Eq. m which are diagonal in 
family space {i.e. with j = k) also vanish since lijk is 
real in this case. 

The equilibrium conditions including non-perturbative 
electroweak reactions, for temperatures below the mass 
scale of superpartners and the critical temperature Tc of 
the electroweak phase transition, yield a relation which 
allows us to find the final baryon asymmetry ns/s in 
terms of the B — SLi asymmetries ns-sLi/s, where nx 
is the density of the quantum number X and s is the en¬ 
tropy density. Note that, in this regime, all three quan¬ 
tum numbers B — ZLi are conserved. This is, however, 
not true for temperatures above the scale of sparticle 
masses, which is taken to lie above the critical temper¬ 
ature of the electroweak transition. We assume that at 
least one of these quantum numbers, which we designate 
as B — 3 L 3 , is conserved at all temperatures after re¬ 
heating (see Sec. IVHII . Then, in a temperature range 
just above the SUSY threshold where B — 3 Li _2 are vio¬ 
lated by processes involving sparticles, we can determine 
ns-sLi j/s in terms of ub-sls/s. Assuming continu¬ 
ity of nB- 3 Li/s as the temperature crosses the SUSY 
threshold, we can express ns/s in terms of nsshs/s 
(see Sec. IV Hill . The result can be written as 


riB _ j^ riB-sLs 
s s 


( 22 ) 


Then if we imagine the inflaton to decay instantaneously 
out of equilibrium creating initial L 3 lepton number den¬ 
sity ns 3 _init, we have 


ns 

s 


= —3k 


n'Ls.iiiit 


S 



9k Treh 

-yes - 

4 TOinf 


(23) 


a = -mf„f Im iijk Y ] ■ 

j,k t,t' 


(19) 


using the standard relation ninf/s = (ns -\- ng)/s = 
3Ti.eh/4niinf for the inflaton number density. The reheat 
temperature is given by 


Here, the indices t, t' can be thought of as indices which 
refer to the topology of the diagrams in Fig. [3 Namely 
t = trg corresponds to the vertex diagram in Fig. |3i, 
while t = self to the self-energy diagram in Fig.l^b. We 
find that 


^trg ^ ^ 

^ IGTT^rnfjjf 




( 20 ) 


where mp ~ 2.44 x 10^® GeV is the reduced Planck scale, 
and g* counts the effective number of relativistic degrees 
of freedom taking account of the spin and statistics and 
is equal to 228.75 for the MSSM spectrum. 


where Li 2 is the dilogarithm [^, and 


VI. EFFECTS OF R-SYMMETRY VIOLATION 


pself _ _± _ ‘ 

^ 47r - M| ’ 


(21) 


We will now explore the possible phenomenological 
and cosmological consequences of the explicit violation 
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of U(l)i{ in the superpotential. Our models contain a 
Z 2 matter parity symmetry under which all the matter 
(quark and lepton) superfields change sign. This sym¬ 
metry is actually the subgroup of U(1 )_r and is left 
unbroken by the soft SUSY-breaking terms. Matter par¬ 
ity combined with the Z 2 fermion parity under which 
all fermions change sign yields R-parity, which, if unbro¬ 
ken, guarantees the stability of the LSP. In the present 
case, however, matter parity is violated along with the 
U(l)jj by the last explicitly displayed term in the RHS 
of Eq. (P) (or Eq. Q), which is necessary for generating 
the observed BAU. Consequently, R-parity is explicitly 
violated and the LSP can decay rapidly, rendering it un¬ 
suitable for dark matter and leading to distinctive collider 
signatures. 

The LSP could decay into a pair of electroweak Higgs 
bosons and a lepton if it contains a Higgsino component. 
The dominant diagrams are constructed from the 11(1)^- 
and R-parity-violating Yukawa vertices with 

the fermionic connected to the fermionic v'j of the 
Yukawa couplings yuij{lih 2 )Vj via a mass insertion. 
For the values of the parameters considered here (see 
Sec.im. we find that the resulting LSP life-time can be 
as low as about 10“^ sec. However, it is easy to block 
kinematically these decay channels of the LSP by taking 
its mass to be smaller than twice the mass of the lightest 
Higgs boson, which is very reasonable. Even then, the 
LSP could decay to a lepton and an electroweak Higgs 
boson with similar life-time. The relevant (one-loop) di¬ 
agrams may be obtained by connecting the two exter¬ 
nal Higgs lines of the previous diagrams to a trilinear 
Higgs vertex. These two-body decay channels can though 
also be blocked kinematically if the LSP mass is taken 
smaller than the lightest Higgs boson mass. In this case, 
the dominant diagrams for the LSP decay may be con¬ 
structed from the above one-loop diagrams by coupling 
the external Higgs line to a pair of fermions (excluding 
the t-quark). So, we obtain one-loop diagrams leading 
to three-body decay of the LSP. They involve an extra 
small Yukawa coupling constant which together with the 
one-loop factor yields a suppressed decay rate. The life¬ 
time can be of order 10® sec or higher and can be further 
enhanced if the LSP is predominantly a gaugino. Thus, it 
may be possible to rescue the LSP as dark matter candi¬ 
date, especially if its Higgsino component is suppressed. 
Alternatively, if the LSP decays fast, we would have to 
ensure that its life-time were less than 1 sec to avoid con¬ 
flict with primordial nucleosynthesis (see e.g. Ref. j^L 

Besides the LSP decay, we could also have other low 
energy processes which violate R-parity. The superpoten¬ 
tial terms which break R-parity involve at least one su¬ 
perheavy field. On integrating out the superheavy fields, 
we generically obtain effective R-parity-violating opera¬ 
tors involving only MSSM fields. If these operators have 
dimension five or higher, they do not lead to detectable 
processes since they are suppressed by some powers of 
Mj. On the contrary, operators with dimension four such 
as the effective scalar vertex hih 2 h 2 l* which originates 


from the superpotential couplings CiC'f(hih 2 ), yi(lih 2 )v'[ 
can lead to low-energy R-parity-violating processes which 
may be detectable in the future colliders (see Sec. HU). 

As in any leptogenesis scenario with RHNs, we must 
ensure that the primordial Li asymmetry (for at least 
one value of i) is not erased by lepton-number-violating 
scattering processes such as klj —> /i 2^2 or lih 2 h^l*^ 
at all temperatures between T^eh and about 100 GeV. In 
our case, due to the presence of the R-parity-violating 
superpotential terms, there exist some extra processes of 
this type such as hilj /i 2^2 or / 11/12 ^ h^l*, which 
come from diagrams similar to the ones mentioned above 
for the fast LSP decay. In addition to all these processes 
which correspond to effective operators of dimension five 
(or higher), we also have dimension-four processes which 
violate R-parity (and lepton number) such as the process 
hih 2 h 2 li derived from the effective four-scalar vertex 
in the previous paragraph. We will return to this issue 
in the next section. 

As already mentioned, the classical flatness of the infla¬ 
tionary trajectory in the limit of global SUSY is ensured, 
in our case, by a continuous R-symmetry enforcing a lin¬ 
ear dependence of the superpotential on S. This is re¬ 
tained lU even after supergravity corrections, given a 
reasonable assumption about the Kahler potential. The 
solution m to the y, problem is also reliant on the R- 
symmetry. These properties are not affected by the ex¬ 
plicit R-symmetry breaking considered here. Note that, 
in our scheme, some R-symmetry-violating couplings are 
included in the superpotential and some not. Thanks to 
the non-renormalization property of SUSY, this situation 
is stable under radiative corrections, but it may be con¬ 
sidered unnatural since there is no symmetry to forbid 
the missing terms. 


VII. PRESERVATION OF R - 3L3 

In this section, we will discuss the conditions for at 
least one quantum number B — iLi to be preserved from 
reheating all the way down to the electroweak phase tran¬ 
sition. This is important since, in the opposite case, the 
primordial lepton asymmetry will be erased. Note that 
the non-perturbative electroweak sphaleron effects pre¬ 
serve all these three quantum numbers. So, violation 
of some B — 3Li can only take place if the correspond¬ 
ing Li is violated perturbatively. (Recall that B is con¬ 
served to all orders in perturbation theory and the non- 
perturbative QCD instanton interactions conserve B and 
all the Li’s.) 

Dangerous Li-violating reactions that could wash out 
any primordial asymmetry in H — 3Li include dimension- 
four scalar interactions arising from the F-term lEycp, 
which violate Li either through the R-parity-violating 
operator C,kV%{hih 2 ) or if they involve two Yukawa cou¬ 
plings yvik{kh 2 )vl and yyjkiljh 2 )vl with i ^ j. Note 
that the latter possibility exists if yi, is not simultane- 
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ously diagonalizable with (the 3x3 matrix with el¬ 
ements yeij)- In the case where the lightest RHN [vf) 
has a mass less than Tj-eh, direct thermal production of 
also takes place. We then have an extra source of 
Li violation from the Yukawa coupling of to li. Fi¬ 
nally, we must consider L^-violating effective dimension- 
five operators arising from fermionic exchange with 
mass insertion (and also without mass insertion in the 
case Treh > Ml). 


A. Dimension-four operators 

As mentioned earlier, individual lepton numbers Li for 
the MSSM superfields can only be defined in the “hatted” 
basis where the charged lepton Yukawa coupling constant 
matrix ye is diagonal. In the SU(2)/j-symmetric case, 
yi, can be diagonalized simultaneously with ye- In the 
non-SU(2)/j-symmetric case, however, this is in general 
not possible. It is though possible to strongly restrict 
the form of fji, by using the remaining freedom of per¬ 
forming unitary rotations in the flavor space of the RHN 
superfields and applying the condition that Li (for some 
particular value of i) is conserved by dimension-four op¬ 
erators. We start by considering the four-scalar operator 
Iih 2 l*h 2 with coupling constant y^ikVljk which, for j ^ i, 
violates Li. In order to retain the Li asymmetry, we must 
impose the condition 

(25) 

k 

for j ^ i. For definiteness, let us choose to preserve the 
L 3 asymmetry, thus i = 3. An appropriate rotation in 
the space can bring the 3-vector y^^k on the 1/3 axis, 
i.e. can make yi/ 31 , yi /32 to vanish. Equation then 
implies that the 3-vectors y^ik and yu 2 k lie in the — V 2 
subspace, i.e. = y^ 2 'i = 0. Using the remaining 

freedom of rotations in this subspace, we can bring yi, 2 fe 
on the V 2 axis, i.e. make ^ 1^21 = 0. So the only non¬ 
vanishing jjvikS are the diagonal ones and yvi 2 - The 
diagonal elements jji of yv can be further made real and 
positive by absorbing their complex phases into the RHN 
superfields, while yvi 2 remains arbitrary and complex. 
For simplicity, we have chosen yyi 2 = 0 so that yi, is 
simultaneously diagonalizable with ye. We still have to 
consider the four-scalar operator /iift, 2 li 2^3 with coupling 
constant ^ 3^3 which also violates L 3 . Therefore, in order 
to retain the L 3 asymmetry, we must further impose the 
condition 

6 = 0 . (26) 

Note that it is by no means necessary for this equation, or 
the other constraints on L 3 -violating couplings, to be sat¬ 
isfied exactly. The dangerous dimension-four operators 
will not come into thermal equilibrium if their effective 
coupling constants are less than about 10“^. However, 
we take them to vanish exactly in order to simplify the 
discussion. 


B. Dimension-five operators 


In Sec. IVHII we will assume that, for a temperature 
range just above the scale of the superpartner masses, 
the Li and L 2 quantum numbers are violated due to 
the four-scalar operators involving C couplings, whereas 
L 3 is perturbatively conserved for all temperatures after 
reheating. The latter is guaranteed for dimension-four 
operators by Eq. but needs to be checked for ef¬ 
fective dimension-five operators with exchange, since 
the RHN mass terms explicitly violate the lepton num¬ 
bers Li. We consider these operators first for the case 
without SU(2)i{ symmetry where the RHN masses are 
all at or above the inflaton mass scale (see Sec. HU), and 
then for the case of SU(2)/{ symmetry. In both cases, 
we estimate the rate of reaction mediated by a chirality- 
conserving RHN exchange as 


F 




-^kT 

^ max(M^, r^) 


(27) 


where Ak is the product of the dimensionless coupling 
constants which corresponds to the diagram with an ex¬ 
change of a RHN with mass Mk. In the case of a chirality¬ 
changing transition with mass insertion on the prop¬ 
agator, the relevant formula is 


F 




AkMk 

^ max(M^,T2) 


2 


(28) 


1. Non-S\J{2)ii-symmetric case 

When all RHNs have masses of the order of mjnf or 
greater, we have Mi > T. Equation (UHli then becomes 

F - \A\^T^, (29) 

where A = Afe/Mfc. For scattering processes of the 
type lih 2 —> ljh 2 , which are mediated by a RHN exchange 
with mass insertion, this is 

2 

Aij = —yi{Mj^c )ijyj — (^ 0 ) 

where, as explained above, we take the neutrino Yukawa 
couplings to be diagonal in the charged lepton flavor ba¬ 
sis (“hatted” basis) and also approximate the up-type 
Higgs VEV to z; ~ 174 GeV (large tan/3 limit). The 
light neutrino mass matrix in the charged lepton basis 
is found as 771 = V*ffiV"^ for given values of the MNS 
matrix V and light neutrino masses rrii. Hence, we find 
that Aij has entries of at most 10“^® GeV“^, assuming a 
hierarchical spectrum of light neutrinos . The ratio of 
the rate of the L 3 -violating scattering processes involv¬ 
ing two Yukawa couplings and a chirality-changing RHN 
exchange to the Hubble rate is thus 

F _ mpr|^ 3 jf _ |2 ^ 4 Q ^ ^q 17 Qgv. (31) 

H 5 
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Therefore, even for reheat temperatures which are as high 
as 3 X 10^° GeV, these scattering processes are well out 
of equilibrium after reheating. 

We also have to take into account scattering processes 
of the type ^ 3/12 ^ hih^ with one Yukawa and one Cj 
coupling, j ^ 3, mediated by a chirality-changing RHN 
exchange. The amplitude for these processes is 

E -^ 3 , (32) 

1=1.2 


For given values of the light neutrino parameters, Q and 
yj, the rate may be found. The worst case occurs in 
the small mi limit for which the rate of these scattering 
processes at T = 3 x lO^^ GeV is estimated as 


r 

H 


3.4 


c 

(mf/GeV) 


2 

190 


c 

(mf/GeV) 


2 


(33) 


where mf* are the (diagonal) Dirac neutrino masses and 
we take Ci = 1^2 with |Ci| = |(^ 2 | = |CI (see Sec. IIXII . 
Hence L 3 is safely preserved after reheating if we have 


nD 


GeV 


>1.8|C|, 


.D 


> 


GeV 


14ICI, 


(34) 


which are easily satisfied within the region of parame¬ 
ter space where the model constraints are fulfilled. Note 
that dimension-five operators mediated by a chirality- 
conserving RHN exchange are even more suppressed as 
one can easily deduce by comparing Eqs. ini) and (EHli- 
So, no L 3 violation by dimension-five operators is encoun¬ 
tered in the non-SU(2)fl-symmetric case after reheating. 


2. SIJ{2)u-symmetric case 

In the SU(2)/i-symmetric case, the neutrino Yukawa 
coupling constants in the “hatted” basis are unambigu¬ 
ously determined to be given by the asymptotic relation 
Vv = (1/f cos/3) diag(me, m^, m^). For the purposes of 
calculating thermal effects near the reheat temperature, 
we can use the asymptotic values of Yukawa coupling 
constants since these constants do not run very much be¬ 
tween the GUT scale and Treh- Our best-fit value of tan (3 
is around 55, thus the numerical values of the yi’s are 
given by iji, = diag(0.00014, 0.028,0.66) approximately, 
taking into account the RG evolution of the charged lep¬ 
ton Yukawa coupling constants to the GUT scale. The 
masses are also determined uniquely by the light neutrino 
masses and the MNS matrix. As mentioned, the RHN 
masses must satisfy M 3 < M^/Ms and M 2 > minif 2. 
The lightest RHN mass (Mi) turns out to be very small 
(of order 10^ GeV) in this case (see Sec. lixl . However, for 
our chosen parameter values, these bounds may be sat¬ 
isfied with a similar choice of light neutrino parameters 
to that for the non-SU(2)ij-symmetric case. Hence the 
entries in the matrix A defined in Eq. are of similar 


size. Thus, for temperatures that are smaller than all the 
RHN masses, the L 3 -violating dimension-five operators 
involving two Yukawa couplings and a chirality-changing 
RHN exchange are out of equilibrium. 

However, we also have to take into account the fact 
that we have a range of temperatures which are higher 
than the lightest RHN mass eigenstate Mi. In this case, 
the contribution of vf to scattering amplitudes will be 
different from its contribution in the previous case. Con¬ 
sidering the propagator for vf exchange with mass in¬ 
sertion, instead of the factor 1/Mi assumed in Eq. 13011 . 
we obtain Mi/T^. But this will evidently give a smaller 
rate than the estimate in Eg. 12911 (which is incorrect in 
this case) if there are no cancellations between contribu¬ 
tions of different RHNs (cancellations appear very un¬ 
likely given the large hierarchy in the masses, and 
this can be checked numerically). Thus, the L 3 -violating 
scattering rate, in this case, is always equal to or smaller 
than the estimate in Eq. (EB, which is still smaller than 
H. Consequently, the Ls-violating scattering processes 
from dimension-five operators which involve two Yukawa 
couplings and a chirality-changing RHN exchange are ex¬ 
pected to be out of equilibrium at all temperatures after 
reheating in the SU(2)/{-symmetric case too. 

We also have scattering through v'^ exchange without 
mass insertion involving two Yukawa couplings. This 
does not violate the total lepton number, but could con¬ 
vert Ls to Lj (j ^ 3). The propagator now varies as 
1/T for T > Ml or T/Ml for T < Mi, which is, in both 
cases, suppressed relative to the dependence of 1/Mi in 
Eq. 13011 . The situation is similar in the case of 3 ex¬ 
change where the discussion at the end of Sec. IVII Ull 
applies. Hence the scattering rate is suppressed relative 
to the rate in Eq. EB, and is again smaller than the 
Hubble rate for temperatures below 3 x 10 ^° QeV. One 
can further show that this conclusion remains valid even 
if we replace one of the Yukawa couplings by a C coupling 
in all cases (with and without mass insertion). 

Since we are considering temperatures above Mi, there 
will be a thermal density of particles and L 3 may be 
violated, say, by direct scattering of and /12 into (see 
also Ref. [^), where the RHN will subsequently decay 
mainly to Zi/i 2 (or lih^)- The rate for this is similar to the 
rate for any dimension-five process through exchange 
without mass insertion (an equivalent process if is put 
on mass shell), except that the coupling constant is sim¬ 
ply ysU^i- The relevant bound is now 

|y 3 U 3 “i P X 4.9 X 10^^ GeV < T, (35) 

which, for T 10^ GeV, implies 

|C/|i|<10-^ (36) 

This bound may in fact be respected for some parameter 
choices, but we will see that the SU(2)/{-symmetric case 
is ruled out by other considerations, namely the smallness 
of off-diagonal elements of 17“. 
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VIII. CHEMICAL POTENTIALS AND ns/s 

In order to find the final baryon asymmetry, we con¬ 
sider the evolution of the net number densities of all the 
relativistic species between the time of reheating and the 
time at which all _B-violating interactions are out of equi¬ 
librium and the baryon number is effectively conserved. 
At any given temperature, we determine which reactions 
(including the non-perturbative ones, i.e. the sphaleron 
or QCD instanton effective vertices) are in equilibrium. 
This is the case if the rate of reaction normalized to the 
particle number density is greater than the Hubble rate. 
Then one can deduce which additive quantum numbers 
are conserved by these reactions. The equilibrium num¬ 
ber density n of a relativistic species minus the number 
density h of its antiparticle is given, in terms of its chem¬ 
ical potential by the standard formula 

n — n _ I5g j 2 (bosons) 
s 1 1 (fermions) 

where g is the number of internal degrees of freedom of 
the species. For each reaction in equilibrium, the sum of 
chemical potentials of the reaction products is equal to 
the sum of chemical potentials of the reactants. Using 
the equilibrium equations, one can express 1 ^ the 
net number density of any species as well as the density 
of any quantum number in terms of a small set of chemi¬ 
cal potentials. Eliminating these chemical potentials, one 
can determine the densities of all quantum numbers given 
the densities of the conserved ones. When some reactions 
go out of equilibrium, one may find the new densities of 
the quantum numbers which cease to be violated by im¬ 
posing that the densities remain unchanged. We use the 
convention that the same symbol denotes the species 
and its chemical potential. For fermions, we count each 
left-handed species separately, hence a non-zero chemi¬ 
cal potential for the gluino g (say) implies an asymmetry 
between left- and right-handed gluinos. 

We will first perform the chemical potential analysis 
in the period which follows reheating and for temper¬ 
atures higher than the mass scale of the superpartners 
and the critical temperature Tc for the electroweak phase 
transition. In this regime, the SM gauge group is un¬ 
broken and we will assume that the spectrum coincides 
with the MSSM one. This is certainly the case in our 
model provided that all the RHN masses are larger than 
Treh- In the SU(2)/{-symmetric case, however, this anal¬ 
ysis will only hold for temperatures which are lower than 
the lightest RHN mass, which is smaller than Treh- Gauge 
interactions imply that the gluons, the photon and 
the B gauge boson have vanishing chemical potential at 
all times, as do the bosons above the electroweak 
phase transition. The interactions of squarks with all 
three types of SM gauge bosons and the corresponding 
gauginos also ensure [33 that all gauginos have the same 
chemical potential, which we denote by g. 

Yukawa interactions originating from the superpoten¬ 



tial are found to be in equilibrium above the electroweak 
scale, implying the relations 

= —q — h 2 , 

= —q — hi, 

el = -k-hi. (38) 


Due to the Cabbibo-Kobayashi-Maskawa mixing, the 
quark chemical potentials are independent of the gen¬ 
eration number. In the leptonic sector, the situation is 
different: due to the very small scale of masses, neutrino 
oscillations are not |29| in equilibrium at the tempera¬ 
tures of interest; and, as has been shown, interactions 
mediated by massive RHN exchange may not be either. 
Thus the charged lepton and neutrino chemical potentials 
may differ between the three generations in the “hatted” 
basis, where the lepton Yukawa couplings are diagonal. 
Note that, although the chemical potentials k (or a) 
and are, throughout, in the “hatted” basis, we sup¬ 
press the hats for simplicity of notation. 

SUSY-breaking soft terms mediate transitions which 
are in equilibrium for temperatures below about 
10^ GeV and above the scale of soft scalar masses (for 
simplicity, we will assume a common SUSY threshold). 
In this regime, the gaugino chemical potential vanishes 
and all members of a supermultiplet have equal chemical 
potential. At higher temperatures, we have the follow¬ 
ing relations between the components of the chiral 
superfields: 


q = q + 9, 

= u'^ + g, 
d‘^ = (f+g, 

g -i 

= e’^ + g, 
hi,2 = hi,2-g. 


Then the baryon 
ns _ 
s 

riLj _ 
s 


and lepton number 


15 

15 

An^g^T 


(4? + 

[3(3^i -f hi 


densities are 
hi + ^ 2 ), 

) + , 


(39) 


(40) 


where Ng is the number of generations. The condition 
that the electric charge density should vanish is 


? = 3[2iVgg-2^Z, + (2iVg + l)(h2-hi)] = 0. 

(41) 

We must now also consider the SU(2)i and SU(3)c 
non-perturbative interactions (’t Hooft vertices), which 
involve all the left-handed fermions transforming non- 
trivially under the corresponding group and are unsup¬ 
pressed at these temperatures. They give the following 
relations: 


SU(2) i, : 3Ngq ^ ^ h hi h2 2^ — 0, 

i 

SU(3)e : Ng{hi + h 2 )- 6 g = 0 , (42) 







13 


which together with Eq. m yield 

J_ 7 _ 2(-/Vg + 3) 
3iV„ ^ * 


q = - 


hi = - 


ho = 


3iV2 


3(27Vg + 1 ) 
4 




leiVg + 3 
3iVg(2iVg + 1) 


3(2iV, + 1) 




2QNg + 15 

iNgilNg + 1 ) 


ff- (43) 


So, we are left only with four independent chemical po¬ 
tentials (Zi, g), in terms of which the baryon and lepton 
number densities are given by 


riB 

s 

riLj 

s 


15 


47r2g*r 




9k- 


lONg - 24 


2 Ng + l^ 
^ 3 


■T,‘i 


ANl + img + 3 
ZVs(21Vg + l) 


g 


(44) 


which hold regardless of the details of the model. 

The scalar potential term \dW/di'^\^ gives rise to op¬ 
erators mediating the transition hiho <-*■ hok for i = 1 , 2 
(recall that we imposed the condition that Ca = 0)- These 
reactions are in equilibrium just above the superpartner 
threshold and imply that = —g + hi. Moreover, as 
already mentioned, below a temperature 10^ GeV and 
above the superpartner mass threshold, g vanishes. So, 
just above the SUSY threshold, we have only one inde¬ 
pendent chemical potential (say I 3 ) and Eq. (IHIl yields 


riB 

s 

U-Li 2 
S 

nL3 

s 


15 -12(2iVg-M) 

in^g^T 6Ng + ll 

15 -48 , 


3) 


47r2g*r 6Ng + ll 

15 3(18A^g-f29) 

in^g^T 6Ng + 11 


(45) 


Eliminating Z 3 , we can express all the quantum numbers 
in terms of one of them. In particular, we hnd that 


riB 

s 

«B-3Li,2 

S 

which also gives 

nB-3Li 


Y 


4(2iVg -I- 1) ub-sLs 
6Ng + 11 s’ 
4(21Vg - 11) nB-3L, 
62Ng -h 91 s 


3(26fVg -I- 1) nB- 3 L 3 
62Ng -h 91 s 


(46) 

(47) 

(48) 


Note that B — ZL^ is the only conserved quantum number 
just above the SUSY threshold. It is actually the only 
quantum number which remains conserved at all times 


after reheating and, thus, the corresponding asymmetry 
is equal to —3 times s. Equation 14(111 then yields 

nB _ 12(2Ng + 1) n_L 3 

62Ng -h 91 s ■ ^ ’ 

This gives the correct result for the final baryon number 
density provided that i) the electroweak phase transi¬ 
tion is strongly first order so that no S-violating SU(2 )l 
sphaleron interactions are in equilibrium in the broken 
phase, and ii) the scale of the superpartner masses is 
smaller than the critical temperature Tc, i.e. the final 
value of ub/s is fixed at a temperature at which reac¬ 
tions involving superpartners are in equilibrium. 

However, if superpartners decouple while sphalerons 
are still active, one must consider the equilibrium equa¬ 
tions without the superpartner contributions, and per¬ 
form a non-SUSY analysis [^. In this case, we have 
two possibilities: either the electroweak phase transition 
is strongly first order (i.e. the sphalerons freeze out very 
quickly after the transition) and the final baryon number 
is fixed in the unbroken phase, or it is second order (or 
weakly first order) and the sphalerons continue to be in 
equilibrium in the broken phase. Recent data indicate 
that the weakly first-order transition is more likely (see 
e.g. Ref. [s^)- We will thus analyze this case in detail. 

With superpartners and heavy Higgs bosons decou¬ 
pled, we are left with the following chemical poten¬ 
tials which a priori may be non-zero: u, d (correspond¬ 
ing, respectively, to the up- and down-type component 
of the left-handed SU(2 )l doublet q), Ci, e°, vt, 

h (corresponding to the lightest Higgs boson), and W~. 
For T > Tc, the chemical potential of the W~ gauge bo¬ 
son vanishes, while below the W~ boson may have 
a non-zero chemical potential, but h = 0 |3^ . This 

immediately fixes the SU(2)j;, singlet chemical potentials 
as u'^ = —u, = —d, el = —e^. Gauge interactions im¬ 
ply d = u + W~, et = + W~, leaving us, for T < Tc, 
with Ng + 2 independent chemical potentials u, Vi, W~. 
The condition that the electric charge vanish leads to the 
relation 

NgU-Y,., - {2Ng + 3)IU- = 0, (50) 


which determines W~ in terms of u and i^i. The baryon 
and lepton number densities can be easily calculated and, 
after eliminating W~ by using Eq. CT . are given by 


riB 

s 

riLj 

s 


15 


2N„ 


47r2g*r 2iVg-h 3 


{5Ng + 6 )u-Y> 


15 

47r2g*r 


^3;zi + 2 


- Ej \ 
2Ng + 3 )' 


(51) 


Now consider the SU(2)i sphaleron interaction. Its equi¬ 
librium condition implies that Ng(u + 2d) + Ei should 
vanish, which fixes u in terms of Vi. 


u = — 


3E^^» 

fVg(8fVg+9)- 


(52) 
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Equation itCT) then yields 

UB 15 4(2Afg + 3)y^ 

s ~ 8Ng + 9 

ns-ZL. = 15 , 4(2jV, - 3) 

s 47rV7’ \ 8 iVg + 9 



(53) 


The three quantum numbers B — ZLi are separately con¬ 
served after superpartners decouple, because reactions 
violating Li ^2 are out of equilibrium. We solve for Vi 
in terms of ns-sLi and substitute back into the baryon 
number density to obtain 

n-B _ 4(2A^g -I- 3) 

s “ 3(321Vg -b 15) ^ s ■ ^ 

From continuity, the constant values of the asymmetries 
n-B-SLi/s for temperatures lower than the SUSY thresh¬ 
old should be identical to the values of these asymmetries 
just above this threshold where the superpartners are still 
in equilibrium and 5 = 0. In this regime, B — 3 Li .2 are 
violated and the values of the corresponding asymmetries 
are expressed in terms of the conserved B — 3 L 3 asymme¬ 
try. Assuming that the SUSY threshold lies higher than 
Tc, the B — 3Li 2 asymmetries are given by Eq. 63 . 
Then, applying Eqs. 6HI> and we find that the final 
baryon number density is given by 


data 0 , *.e. — 0.1126, and putting tolsp = 

100 GeV, Eq. 63 yields T^eh 3.04 x lO^o GeV. This 
high value of Tj-eh, chosen here to maximize the BAU, 
can well be 0 compatible with the gravitino constraint 
iH provided that the branching ratio of the gravitino 
decay to a photon and a photino is adequately small. 

The R-symmetry-violating effective coupling constants 
Cl, C 2 are taken to satisfy the relation Ci = *-e. the 

CP-violating relative phase between them is taken equal 
to 7 r/ 2 . The magnitude of these coupling constants is 
maximized, i.e. |Ci| = IC 2 I = M/Mg- These choices 
maximize the generated baryon asymmetry. In order to 
preserve perturbativity and to have a good effective field 
theory, we require M < O.IM 5 . 

As already mentioned, the light neutrino mass-squared 
differences and mixing angles are generally taken to lie 
in the 2cr confidence intervals in Ref. 1^. In the non- 
SU(2)/j-symmetric case, in particular, we take them to 
coincide with their best-fit values in Ref. 6 m\i = 

8.1 X 10-’5 eV^ = 2.2 x IQ-^ eV^ sin^0i = 0.5, 
sin^ 6*2 = 0, and sin^ 6 ^ = 0.3. For simplicity, in the MNS 
matrix V, we set the Dirac phase S and the Majorana 
phases a and j3 to zero. This highlights the fact that 
our mechanism works in the absence of any CP-violating 
phases in the standard (R-parity-conserving) couplings. 


A. Non-SU( 2 )ii-symmetric case 


ns _ 4:{26Ng + l){2Ng + 3) bb-sl^ 

~ ~ (62Ag-b91)(32Ag-bl5) s 

_ 12{26Ng + l){2Ng + 3) ul, ,init 

~(62Arg-b91)(32Ag-bl5) s 

2844 init /'reel 

~10249 s ^ ^ 

for Ng = 3. 

IX. NUMERICAL RESULTS 

We make the assumption that the LSPs are long-lived 
and their thermal production is negligible, which holds 
in many cases. So, the LSPs which possibly contribute 
to the cold dark matter (CDM) in the present universe 
can originate solely from the late decay of the gravitinos 
which were thermally produced at reheating. Their relic 
abundance is 0 

where Dlsp is the relic density of the LSPs in units of the 
critical density of the universe, h is the present Hubble 
parameter in units of 100 km sec“^ Mpc“^, and ttilsp 
is the LSP mass. We further assume that the CDM in 
the universe consists exclusively of LSPs. Taking the 
best-fit value for the CDM abundance from the WMAP 


We will first examine the non-SU(2)/j-symmetric case. 
We fix the parameter k to the value 10“'*. The cos¬ 
mic microwave background explorer (COBE) value of 
the quadrupole anisotropy of the CMBR {{6T/T)q ~ 
6.6 X 10“®) [0 is then reproduced for A ~ 3.48 x 10“"* 
(> K as it should) and M ~ 5.63 x 10*® GeV. Thus, 
TOinf — 7.96 X 10** GeV. Note that the values of A and 
TOinf are consistent with the value of T^eh (see Eq. 63 )■ 
The spectral index of density perturbations comes out 
practically equal to unity. 

For definiteness, we set tan/3 = 50, although in this 
case it scarcely affects the calculation of baryon asym¬ 
metry. We are still free to vary Ms, mi and the three 
neutrino Dirac masses m^. As previously mentioned, 
we require that all three RHN masses be greater than 
TOinf /2 and less than M'^/Ms, i.e. Mi > minf /2 and 
M^Ms < M^. By placing one RHN mass Mj near 
the inflaton pole Mj = minf, the contribution of the 
corresponding self-energy diagram to the baryon asym¬ 
metry is resonantly enhanced (see Eq. 12111 1. However, 
we regard this as an unnatural fine-tuning since it ap¬ 
parently cannot be ensured by any symmetry, and re¬ 
ject solutions where any Mj is less than 10 % from the 
pole, i.e. \Mj — minf|/winf < 0.1. We also impose the 
requirement that L 3 is not violated by any dimension- 
five processes and that the observed baryon asymmetry 
ub/s ~ 8.66 X 10“**, derived from the recent WMAP 
data 01, is reproduced. Under all these restrictions and 
for fixed values of mf, we find the value of mi which 
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FIG. 3: The contours with fixed maximal Ms in the mf* — plane for = 30 GeV. The cyan, green, magenta, red and 
blue contours correspond, respectively, to Ms ~ 6 , 7, 8 , 9 and 10 in units of 10^® GeV. In black and white, contours with 
higher Ms appear to be darker. The red/dark shaded area is excluded by both the requirements that, for the maximal Ms, 
Ml > and L 3 is conserved by dimension-five processes, while the blue/lightly shaded one only by the latter. 


maximizes Mg (and thus minimizes ICil and |C 2 |)- We 
have chosen to maximize Mg in order to obtain values as 
close as possible to the standard scale (f« 5 x 10^^ GeV) of 
the weakly-coupled heterotic string. We keep only values 
of maximal Mg which exceed lOM ~ 5.63 x 10^® GeV so 
that the perturbativity requirement is satisfied. 

In Figs. 13 and ^ we present the contours with fixed 
value of this maximal Mg and the corresponding mi re¬ 
spectively in the mf* — plane for = 30 GeV. In 
Fig.El we include only contours with fixed maximal Mg 
which satisfies the perturbativity limit. However, this 
limit is not necessarily satisfied everywhere on the con¬ 
tours with fixed mi which we show in Fig. 0] It holds 
only on the parts of these contours which lie in the area 
covered by the contours in Fig. 13 We observe that the 
contours with fixed maximal Mg show the existence of 
two main almost horizontal mountain ranges not far from 
each other and a lower almost vertical mountain range at 
low values of mf*. The contours with fixed mi, on the 


other hand, form a single horizontal mountain range and 
a secondary almost vertical mountain range at higher val¬ 
ues of m 2 ■ At low values of (and mf*), we find that 
Ml < for a wide range of choices of mi and our 

scenario is irrelevant since the inflaton decays directly 
into RHNs. Also, the ratios Q/Vj (j = 2) acquire large 

values, which result in the L 3 asymmetry being washed 
out by dimension-five processes involving a C coupling. 
The red/dark shaded areas in Figs. 0 and 0 are excluded 
because, for the maximal value of Mg, Mi < mini 12 and 
L 3 is violated; the blue/lightly shaded areas, on the other 
hand, are excluded only because of L 3 violation. The 
value of the maximal Mg in these areas is much smaller 
than the perturbativity bound, thus these areas are not 
considered acceptable in any case. We have also con¬ 
structed the contours in the mf* — mf* plane with fixed 
Mg or mi for m^ = 10 GeV and mf* = 20 GeV. We 
find that they reveal a very similar structure. 

In order to understand the structure of these con- 
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FIG. 4: The contours with fixed mi corresponding to the maximal Ms in the mf* — plane for = 30 GeV. The cyan, 
green, magenta, red and blue contours correspond, respectively, to m\ = 3, 5, 10, 50 and 100 in units of 10“"^ eV. In black and 
white, contours with higher mi appear to be darker. The red/dark shaded and blue/lightly shaded areas are as in Fig. 



tour plots, let us consider some specific values for the 
three Dirac neutrino masses mf'. For instance, let us fix 
mf = 7 GeV, = 30 GeV and take five characteristic 
choices of (4, 4.35, 4.6, 4.95, 5.3 GeV). Then com¬ 
pute the value of Ms which is required in each case to 
obtain the correct ns/s as a function of mi. The result 
is depicted in the top left panels of Figs. I5I9I bv a solid 
line. Moreover, we find that the constraint M 3 M 5 < 
is satisfied to the right of the dashed line in these panels. 
In the other three panels in each of these figures, we show 
the corresponding values of Mi, i = 1,2, 3, as functions of 
mi- The shaded bands are excluded by the requirement 
that |Mi — minfl/minf > 0 . 1 . 

To appreciate the meaning of these figures, we must 
first observe that, as our numerical findings show, the 
main contribution to the baryon asymmetry comes from 
the interference of vertex with self-energy diagrams, 
while the interference between self-energy diagrams is 
subdominant. The heaviest RHN mass M 3 is much larger 
than Ml, M 2 and m^f as one can see from Figs. I5I9I 


So, the contribution to the baryon asymmetry from dia¬ 
grams with 1/3 exchange is suppressed. Moreover, in the 
interesting range of parameters where the maximal Ms 
is achieved. Mi is nearer the inflaton pole than M 2 is. 
As a consequence, the dominant self-energy diagram is 
the one with vf exchange. It also turns out from our nu¬ 
merical study that the dominant vertex diagram in the 
interesting range of parameters is the one with i >2 ex¬ 
change. So, the baryon asymmetry is dominated by the 
interference of the self-energy diagram with vf exchange 
with the vertex diagram with 1/2 exchange. 

From Figs. EEl we see that, for toi ^ 0 , the masses 
Ml and M 2 acquire constant asymptotic values, while 
M 3 keeps increasing. In consequence, the kinematic fac¬ 
tor in the dominant contribution to ns Is, which origi¬ 
nates from the “stripped” diagrams, is constant in this 
limit. The prefactor originating from the product of cou¬ 
pling constants contains [/31 and 1/32 which connect the 
i >3 state with the 2 mass eigenstates. We find that, as 
mi decreases, these elements of first rise sharply and 
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FIG. 5: The value of Ms required to obtain the correct baryon 
asymmetry as a function of mi (solid line in the top left panel) 
for mf = 7 GeV, mf = 4 GeV and mf = 30 GeV. The 
constraint M 3 Ms < M^ holds to the right of the dashed line 
in this panel. The other three panels show the masses Mi 
(i = 1,2,3) as functions of mi for the same values of the 
Dirac neutrino masses. Regions excluded by the requirement 
that |Mi — minf l/minf > 0.1 are shaded. 




FIG. 6: As in Fig. El but for mj* = 4.35 GeV. 


then approach their constant asymptotic values. Their 
rapid increase is compensated by an appropriate increase 
of Ms, which enters the prefactor through |^i| and IC 2 I, so 
that the baryon asymmetry remains equal to its WMAP 
value. This explains the fact that, as mi decreases, 
Ms generally exhibits a fast increase before entering its 
asymptotic plateau (see solid line in the top left panels of 
Fig'S. I5I0II . This phenomenon generally helps us to obtain 
larger maximal values of Ms. 

From these figures, we also see that M 2 and M 3 (see 
two bottom panels) remain effectively unaltered as we 
vary for fixed mf* and m^. This explains also the 
fact that the constraint M^Ms < M^, as depicted in the 
top left panels (dashed line), is practically unaffected by 


FIG. 7: As in Fig. El but for mf = 4.6 GeV. 


changes in the value of mf. The Mi though is enhanced 
as mf increases. For mf = 4 GeV, we see from Fig. El 
that the asymptotic value of Mi (as mi ^ 0 ) lies well 
below the inflaton pole. So, there is no enhancement of 
the dominant “stripped” self-energy diagram produced 
by this pole. As a consequence, the values of M 5 yielding 
the correct baryon asymmetry are relatively low. So, 
the maximal Ms, which is achieved at M 3 Ms = M^ 
in this case, is also quite low and corresponds to a low 
value of mi (^ 10“"* eV). This explains the fact that 
this particular case lies well below the lower horizontal 
mountain range in Fig. El and the horizontal mountain 
range in Fig. El 

Let us now discuss how the situation changes as we in¬ 
crease mf with mf and mf fixed. For mf = 4.35 GeV, 
Fig.Elshows that the asymptotic Mi, which is larger than 
in the previous case, gets closer to the pole and thus the 
dominant “stripped” self-energy diagram is enhanced, 
leading to a larger maximal Ms- This is again achieved 
at M^Ms = M^ and corresponds to a somewhat bigger, 
but still quite small value of mi. This case lies near the 
brow of the lower mountain range in Fig. El but still be¬ 
low the mountain range in Fig. El Increasing mf further 
to the value 4.6 GeV, we see from Fig.Elthat the asymp¬ 
totic Ml gets very close to the pole and, thus, Ms is fur¬ 
ther enhanced. However, Mi now enters into the band 
excluded by the requirement |Mi — minf|/minf > 0 . 1 , 
which in the top left panel extends to the right of the 
line M 3 M 5 = M^. So, in this case, the maximal Ms 
is achieved at the boundary of this excluded band and 
not at M^Ms = M^. It is thus quite low and corre¬ 
sponds to a considerably larger mi. This case lies in the 
valley extending between the two horizontal mountain 
ranges in Fig. El and very near the brow of the horizontal 
mountain range in Fig. El In Fig. El we present the case 
mf = 4.95 GeV. We see that the asymptotic Mi now lies 
above the excluded band, but still close to the pole. The 
whole excluded band in the top left panel now lies to the 
right of the line M 3 M 5 = M^, and thus the maximal Ms 
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FIG. 8: As in Fig. but for = 4.95 GeV. 


is again achieved on this line and corresponds to a low 
mi. Therefore, this case lies above the brow of the hori¬ 
zontal mountain range in Fig.^ Actually, at some value 
of m^, the point at which the maximal Ms is achieved 
jumps from the right boundary of the excluded band to 
the line M 3 M 5 = M^. So, mi drops suddenly to low val¬ 
ues. This discontinuity explains the absence of the upper 
branches of the contours with mi = 10“^, 5 x 10“^ and 
10“^ eV. The enhancement of the self-energy diagram 
due to the pole yields the sharp peak of Ms within the 
excluded (shaded) band. The tail of this enhancement 
leads to a large value for the maximal Ms- So, this case 
is near the brow of the upper mountain range in Fig. 13 
For = 5.3 GeV, the asymptotic Mi moves away from 
the pole and the excluded (shaded) band in the top left 
panel further to the right (see Fig.0. Consequently, the 
maximal Ms is still achieved on the line M^Ms = M^, 
but drops to lower values. This case is well above the 
upper horizontal mountain range in Fig. |3 and the hori¬ 
zontal mountain range in Fig.^ 

The secondary vertical mountain range in Fig. 13 can 
be understood as follows. As mf decreases for fixed and 
large values of mf, the mass Mi also decreases. As a 
consequence, the band excluded by the requirement that 
I Ml — minf l/minf >0.1 moves to the left in the mi — Ms 
plane approaching the curve M^Ms = M^ and the max¬ 
imal Ms, which is achieved on this curve, increases since 
the inflaton pole gets nearer. At some small value of 
mf, the left boundary of this band touches the point 
corresponding to the maximal Ms, which thus reaches its 
largest possible value for the chosen value of ■ This 
yields the brow of the vertical mountain range in Fig. 0 
The corresponding value of mi is, however, still relatively 
small. For smaller values of mf*, the band moves further 
to the left and the point of maximal Ms jumps suddenly 
to the right boundary of the band. Thus, the maximal 
Ms drops, while the corresponding mi rises sharply form¬ 
ing the vertical mountain range in Fig. 21 This sharp dis¬ 
continuity in mi explains the absence of the right branch 





FIG. 9: As in Fig.^l but for = 5.3 GeV. 

of the contour with mi = 10“^ eV. As the value of mf* 
decreases further, the excluded band shifts further to the 
left and the maximal Ms and the corresponding mi drop 
smoothly. 

Despite the abundance of adjustable parameters in the 
non-SU(2)fl-symmetric case, it turns out not to be easy 
to achieve the observed baryon asymmetry. In partic¬ 
ular, we find that, for central values of the neutrino 
mass-squared differences and mixing angles, the string 
scale Ms is restricted to be lower than about 10^^ GeV. 
So, the standard weakly-coupled heterotic string scale 
(« 5 X 10 ^^ GeV) cannot be achieved in this case. How¬ 
ever, lower string scales, such as the ones encountered 
here, can be easily obtained in other string models like 
the stron gly- coupled heterotic string from M-theory (see 
e.g. Ref. |^). Also, we find that it is difficult to gener¬ 
ate an adequate baryon asymmetry with T^eh 10® GeV, 
which is the standard upper bound on T^eh from the grav- 
itino constraint We are thus obliged to allow higher 
reheat temperatures, which is though perfectly pos¬ 
sible provided that the gravitino decay to photons and 
photinos is somewhat suppressed. 

At the upper bound on Ms (« 10^^ GeV), the pa¬ 
rameters |(^i| {i = 1,2) are about 5.63 x 10“® and may 
thus lead to processes observable in the future colliders. 
(For smaller values of Ms, they are even larger.) Indeed, 
as mentioned in Sec. ED the explicit R-parity violation 
in our model, required for leptogenesis, has some low- 
energy signatures coming from dimension-four effective 
scalar vertices. They may typically be the three-body 
slepton decay processes 

h - hM*, (57) 

for i = 1 , 2 , which can easily be kinematically allowed. 
The magnitude of their effective coupling constants (^iy* 
is of order 10 “® near the largest possible value of Ms, 
which is achieved at mf* « 3 GeV and mf « 5 GeV 
for mf = 30 GeV (see Fig. 2)). The corresponding decay 
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rates are then of order 10“® GeV for mass of the decaying 
slepton of order 1 TeV. 

We have seen that the maximal value of Ms, in most 
cases, is achieved at M^Mg = M^, and the correspond¬ 
ing value of mi lies in the range 10“"* — 10“^ eV. How¬ 
ever, for values of the Dirac neutrino masses in the val¬ 
ley between the two horizontal mountain ranges appear¬ 
ing in the Mg contour plots in the mf* — plane for 
fixed (see Fig. 0 ), the maximal Mg is achieved at 
the boundary of the band excluded by the requirement 
I Ml — minf l/winf >0.1 and the corresponding mi is con¬ 
siderably larger reaching values of order 10“^ eV. This 
case, however, yields very low maximal values of M 5 and 
must, therefore, be excluded by the perturbativity re¬ 
quirement. In conclusion, we predict that, for acceptable 
values of Mg, the (in principle observable) smallest neu¬ 
trino mass mi takes values between 10“^ and 10“^ eV. 
Thus our assumption of a hierarchical light neutrino spec¬ 
trum is self-consistent. 


B. SU(2)i{-symmetric case 

In the SU(2)ij-symmetric case, we again put k, = 10“^. 
The COBE value of the quadrupole anisotropy of the 
CMBR is now reproduced for A ~ 3.45 x 10“^ and M ~ 
5.69 X 10^5 GeV, yielding minf 8.05 x 10“ GeV. The 
spectral index is again equal to unity for all practical 
purposes. In this case, tan/3 ~ 55 and the asymptotic 
values of jji are 0.00014, 0.028 and 0.66 for z = 1, 2 and 
3 respectively, as already mentioned. So, in the SU(2)u- 
symmetric case, only Mg and mi remain free. 

To give a chance to this very restrictive case to possibly 
yield acceptable values of the baryon asymmetry, we al¬ 
low all the neutrino oscillation parameters to vary within 
their 2(j confidence interval given in Ref. (for simplic¬ 
ity, we take the Dirac phase (5 = 0). Even then, we find 
that, when the restriction Ca = 0 from preservation of the 
B — 3 L 3 asymmetry is imposed, the resulting values of 
the baryon asymmetry are always some orders of magni¬ 
tude smaller than the observed value for any reasonable 
values of Mg (satisfying the perturbativity requirement) 
and any mi in the allowed range from WMAP data . 

The main reason for this failure is the smallness of the 
off-diagonal elements of C/'^. The existence of a CP asym¬ 
metry requires that vertices be connected to or C ,2 
vertices through internal lines, yielding products of the 
type for fixed values of j and k (j ^ k). 

These products contain at least two off-diagonal elements 
of However, the strongly hierarchical values of iji 
resulting from the SU(2)ij symmetry imply that the off- 
diagonal elements of U‘^ are always small. In addition, the 
requirement from the preservation of the B — 3 L 3 asym¬ 
metry that |17|i| < I0“® further restricts this product of 
matrix entries, such that the overall product of coupling 
constants entering the L 3 asymmetry per inflaton decay 
(£ 3 ) is already less than 10 “®, without considering any 
loop suppression factors. In contrast, in the non-SU(2)ij- 


symmetric case, the values of the off-diagonal elements 
of t/° can be of order 0.1 or even unity. 


X. CONCLUSIONS 

We proposed a scenario of non-thermal leptogenesis 
following SUSY hybrid inflation, in the case where the 
light neutrinos acquire masses exclusively via the stan¬ 
dard seesaw mechanism, i.e. through their coupling to 
heavy RHNs, and the decay of the inflaton to RHN su¬ 
perfields is kinematically blocked or (in the presence of 
a SU(2)/j gauge symmetry) the RHNs with mass smaller 
than half the inflaton mass are too light to generate suffi¬ 
cient baryon asymmetry through their subsequent decay. 
The primordial lepton asymmetry is generated through 
the direct decay of the inflaton into light particles. We 
explored our scenario within the context of two simple 
SUSY GUT models, one without and one with SU(2)i{ 
gauge symmetry, which incorporate the standard version 
of SUSY hybrid inflation. 

The y problem is solved via a U(l) R-symmetry which 
forbids the existence of an explicit /i term, while allows 
a trilinear superpotential coupling of the gauge singlet 
inflaton superfield to the electroweak Higgs superfields. 
After the spontaneous breaking of the GUT gauge sym¬ 
metry, this singlet inflaton acquires a suppressed VEV 
due to the soft SUSY-breaking terms. Its trilinear cou¬ 
pling to the Higgs superfields then yields a /x term of the 
right magnitude. 

The main decay mode of the inflaton is to a pair of 
electroweak Higgs superfields via the same trilinear cou¬ 
pling. The initial lepton asymmetry is created in the 
subdominant decay of the inflaton to a lepton and an 
electroweak Higgs superfield via the interference of one- 
loop diagrams with exchange of different RHNs. The 
existence of these diagrams requires the presence of some 
specific superpotential couplings which explicitly violate 
the U(l) R-symmetry and R-parity. These couplings, 
however, do not affect the exact baryon number conser¬ 
vation in perturbation theory which is implied by the 
R-symmetry. Thus, the only way to generate baryons 
in these models is via a primordial leptogenesis (or via 
electroweak baryogenesis). 

In our analysis, we took into account the constraints 
from neutrino masses and mixing. The requirement 
that the primordial lepton asymmetry not be erased by 
lepton-number-violating processes before the electroweak 
phase transition is a much more stringent constraint on 
the parameters of the theory. We showed that the model 
with SU(2)/{ gauge symmetry is too restrictive to be able 
to generate an adequate baryon asymmetry in accordance 
with these constraints and for natural values of the other 
parameters even if we allow the neutrino oscillation pa¬ 
rameters to vary within their 2cr confidence intervals. It 
is thus ruled out. 

On the contrary, the non-SU(2)/j-symmetric model can 
be viable even with central values of the neutrino mass- 


20 


squared differences and mixing angles. However, we find 
that this model is much more restrictive than the model 
studied in Ref. Ea, which contained SU(2)j;, triplet su¬ 
perfields giving a second contribution to light neutrino 
masses after that of the RHNs. Indeed, in order to gen¬ 
erate the observed BAU, we had to take a larger reheat 
temperature, which is though perfectly acceptable if the 
gravitino decay to photons and photinos is somewhat 
suppressed, and a string scale somewhat smaller than the 
weakly-coupled heterotic string one. Such lower string 
scales are easily obtained in other string models such 
as the strongly-coupled heterotic string from M-theory. 
We also find that the lightest neutrino mass eigenvalue, 
which is an in principle measurable parameter, is re¬ 
stricted to lie in the range 10“"^ —10“^ eV. So, our model 
is consistent with a hierarchical neutrino mass spectrum. 
The explicit breaking of R-parity, which is necessary for 
our baryogenesis mechanism, need not have currently ob¬ 
servable low-energy signatures, although it may have sig¬ 
natures detectable in future colliders. Also, the LSP can 
be made long-lived and, thus, be a possible candidate for 
the CDM in the universe. 

There are many constraints and necessary conditions 
for our mechanism to successfully produce the observed 
BAU (in the non-SU(2)i{-symmetric case). Some of these 
constraints are generic and easily satisfied, but some re¬ 
strict certain a priori adjustable parameters of the model 
to lie within narrow ranges. This is not a technical fine- 
tuning problem, but might be regarded as an undesirable 
feature since we have not presented an underlying theory 
which selects these particular ranges of values. However, 
given the large number of possible fundamental theories, 
it is plausible that parameter values consistent with our 
baryogenesis mechanism will emerge naturally from one 
or many of them. For such underlying theories, which on 
the basis of previous baryogenesis scenarios one might 
think were ruled out, we have shown that in fact they 
may be consistent with data. Furthermore, the scenario 
is especially predictive and is more readily testable com¬ 
pared to baryogenesis models in which parameter values 
are not particularly restricted. 

The conditions for the baryogenesis mechanism to be 
effective involve the following parameters: the dimension¬ 
less (complex) effective coupling constants Q , the lightest 
neutrino mass, the RHN masses, the neutrino Yukawa 
coupling constants (including their complex phases), the 
inflaton mass and the reheat temperature. Recall that, 
for simplicity, the light neutrino mass-squared differences 
and mixing angles are taken to coincide with their best-fit 
values and the complex phases in the MNS matrix are set 
to zero. In addition, the model is restricted to be con¬ 
sistent with cosmological observations determining the 
values of the inflationary parameters and restricting the 
reheat temperature and the mass and couplings of the 
dark matter candidate. 

The first necessary condition is that all RHN masses 
should exceed half the inflaton mass. This is not a fine- 
tuning. Indeed, it may be generic in many classes of 


models. The hierarchical light neutrino mass spectrum 
is similarly restricted to be of the normal (non-inverted) 
hierarchical type. We imposed a condition that RHN 
masses should not be too close (within 10 %) to the infla¬ 
ton mass (Sec. ITxTni specifically in order to avoid regimes 
where the baryon asymmetry is due to a resonant en¬ 
hancement. These regimes would in principle be allowed, 
but the resonance condition appears coincidental since 
the RHN masses and the inflaton mass arise from dif¬ 
ferent sources. We could remove this condition and the 
allowed regions would be slightly larger, but the expla¬ 
nation of the observed baryon asymmetry would be more 
complicated in the resonant regions, hence for simplicity 
we do not consider them. 

The U(l)i{-breaking coupling constants Q are fixed 
within relatively narrow ranges of values. One of them, 
^ 3 , must be quite small (of order 10 “^) or vanish in or¬ 
der not to wash out the produced lepton asymmetry in 
the L 3 direction (Sec. IVIIII . This might, for example, 
be ensured by imposing a symmetry. The other two 
must take values of a few times 10“^. The upper bound 
on them arises from the perturbativity limit on M/Ms, 
while the lower from the requirement to produce suffi¬ 
cient ubIs. We have taken these two coupling constants 
equal in magnitude for simplicity, but they may be un¬ 
equal without much affecting the success of the model. 
Since these values are not particularly tiny and are not 
required to cancel against any other quantity, the restric¬ 
tions represent an in principle testable prediction rather 
than a fine-tuning. The relative complex phase between 

and (C 2 is of order unity. 

The neutrino Yukawa coupling constant matrix (or 
Dirac mass matrix after SU(2)j;, breaking) is taken to 
have small or vanishing off-diagonal elements, when writ¬ 
ten in the “hatted” weak interaction basis where lepton 
family numbers are defined (the charged lepton mass ba¬ 
sis) . Three off-diagonal elements can be set to zero by re¬ 
defining the RHN superfields, two of the remaining three 
must be small or zero to prevent washout of the L 3 asym¬ 
metry (Sec. IVHl . and the remaining one {yui 2 ) is set to 
zero for simplicity. These requirements are not fine-tuned 
since the off-diagonal y^ij ’s may be small as a result of a 
symmetry, as in the quark sector. 

Cosmological parameters that are fixed by cosmo¬ 
logical considerations and data are the number of e- 
foldings of our present horizon scale and the quadrupole 
anisotropy of the CMBR, which together determine M, 
and one of k or A (Sec. The reheat temperature is 
also constrained (Sec. CSI) not to overproduce LSPs from 
gravitino decay, and (less stringently) by the usual grav¬ 
itino constraint on dissociation of light elements. Since 
the lepton asymmetry is proportional to the reheat tem¬ 
perature, we choose the reheat temperature at or near 
its maximum value, which fixes the remaining one of k 
and A. As a consequence, the inflaton mass, which enters 
into the baryon asymmetry formula, is also fixed. These 
values represent neither predictions nor fine-tunings, but 
simply observational constraints on the model. The con- 
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ditions that the LSP be long-lived and the gravitino de¬ 
cays not disrupt primordial nucleosynthesis both restrict 
the spectrum of superpartners, but do not directly influ¬ 
ence the baryogenesis scenario (except through the effect 
of the LSP mass on the maximal reheat temperature). 

The remaining parameters are mf*, the diagonal Dirac 
neutrino masses, corresponding to the neutrino Yukawa 
coupling constants. Although the value of can be 
varied over a range of a few tens of GeV without much 
change, the model is rather sensitive to rn^ 2 i can be 
seen in Figs. 13 and ^ The reason for this is, roughly, that 
the correct primordial lepton asymmetry can be obtained 
only when one RHN mass is not very far from minf. This 
approximate condition results in a surface in the space 
of mf* near which the model is successful (see Sec. El 
for a full discussion). Since mf ’2 protected by chiral 
symmetry, their values are not subject to a technical fine- 
tuning problem, and their rather narrow allowed ranges 
can be said to be a prediction of the model, which is in 
principle testable. It might be argued that it is undesir¬ 


able for the baryon asymmetry to vary steeply when the 
basic parameters of the model are changed by a small 
amount; and in fact we have excluded the “resonance” 
regions where one RHN mass is too close to Winf, where 
this steep variation does occur. Note, though, that we 
do not predict ns/s, but rather use its known value as a 
constraint on the model. 
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